Rashba spin torque in an ultrathin ferromagnetic metal layer 
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In a two-dimensional ferromagnetic metal layer lacking inversion symmetry, the itinerant electrons 
mediate the interaction between the Rashba spin-orbit interaction and the ferromagnetic order 
parameter, leading to a Rashba spin torque exerted on the magnetization. Using Keldysh technique, 
in the presence of both magnetism and a spin-orbit coupling, we derive a spin diffusion equation 
that provides a coherent description to the diffusive spin dynamics. The characteristics of the spin 
torque and its implication on magnetization dynamics are discussed in the limits of large and weak 
spin-orbit coupling. 
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I. INTRODUCTION 

By transferring angular momentum between the elec- 
tronic spin and the orbital, spin-orbit coupling fills the 
need for electrical manipulation of spin degree of freedom. 
Outstanding examples are the electrically generated bulk 
spin polarizationii^ and the well-known spin Hall effect 
(SHE)^"— in a two dimensional electron gas where the 
spin-orbit interaction, particularly of the Rashba-typC)^ 
plays the leading role. Rashba spin-orbit interaction not 
only introduces an effective field perpendicular to the lin- 
ear momentum but also provides the backbone to the 
spin-relaxation through the so-called D'yakonov-Perel 
mechanism^ which is dominant in a two-dimensional 
system. Besides its prominent role in semiconductors, 
Rashba spin-orbit coupling is believed to exist at ferro- 
magnetic/heavy metal as well as ferromagnetic/metal- 
oxide interfaces, in which the inversion symmetry break- 
ing offers a potential gradient empowering the spin-orbit 
coupling. 

Meanwhile, magnetism continuously stimulates the in- 
dustrial and academic appetite. In the pursuit of fast 
magnetization switching, Slonczewski-Berger spin trans- 
fer torque^ employs a polarized spin current instead of 
a cumbersome magnetic field. This celebrated scheme 
demands non-collinear magnetic textures in forms of, for 
example, spin valves or domain wall structures.^ 

In the presence of inversion symmetry breaking (such 
as asymmetric interfaces), a ferromagnetic metal layer as- 
sembles both magnetism and spin-orbit coupling, hence 
offering an alternative switching mechanism i^°'^^ Spin- 
orbit coupling transfers the orbital angular momentum 
carried by an electric current to the electronic spin, thus 
creating an effective magnetic field (Rashba field). As 
long as the effective field is mis- aligned with the magne- 
tization direction, the so-called Rashba torque emerges, 
thus exciting the magnetization. 

Current-driven magnetization dynamics by spin-orbit 
torque has been demonstrated by several experiments 
on metal-oxide based systems J^"— In fact, the Rashba 
torque can be categorized into to a broader family 
of spin-orbit interaction induced torque that has been 
observed in diluted magnetic semiconductors J^"— Re- 



cently, Miron et al.,^^ has demonstrated the current- 
induced magnetization switching using a single ferromag- 
net in Pt/Co/AlOa; trilayers, which further consolidates 
the feasibility of the Rashba torque. The same type of 
spin-orbit coupling induced torque is predicted to im- 
prove current-driven domain wall motionjiiiiS which is 
supported by experimental observations At this stage, 
we are aware of an alternative explanation, as pointed 
out by Liu et al.,^^ in terms of the spin Hall effect (SHE) 
occurring in the underlying heavy metal layer, such as 
Pt or Ta. The distinction between the spin Hall induced 
effect and the Rashba one is discussed in the last section 
of this article. 

In searching for a general form of the Rashba torque 
in ferromagnetic metal layers^ we found an expression 
that consists of two componentsi^An in-plane torque 
(oc m X (y X m)) and an out-of-plane one (oc y x m), 
given y is the in-plane direction transverse to the injected 
current and m is the magnetization direction. Numerical 
solution on a two-dimensional nano-wire with one open 
transport direction has been carried out to appreciate 
the significance of diffusive motion on the spin torque. 
We found that the in-plane component of the torque in- 
creases when narrowing the magnetic wire^^. 

In the present article, we give a full theoretical deriva- 
tion of the coupled diffusive equation for spin dynam- 
ics in a ferromagnetic metal layer and describe the form 
of the Rashba torque in both weak and strong Rashba 
limits. In Sec. [Hi we combine the Keldysh formalism 
and the gradient expansion technique to derive a cou- 
pled diffusion equation for charge and non-equilibrium 
spin densities. To demonstrate that the diffusion equa- 
tion provides a coherent framework to describe the spin 
dynamics, we dedicate Sec. IIIII to the spin diffusion in 
a ferromagnetic metal, which shows an excellent agree- 
ment to early result on the same system. In Sec lIVl we 
illustrate that the absence of magnetism (in our diffusion 
equation) describes the well-know phenomenon of elec- 
trically induced spin polarization. The cases of a weak 
and a strong spin-orbit coupling are discussed in SecfVl 
and Sec. IVIl respectively, where we provide an analytical 
form of the Rashba torque in an infinite medium. In Sec. 
IVIIl we discuss the implication of the Rashba torque on 
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magnetization dynamics as well as its distinction from 
spin Hall effect induced torque. 



II. DIFFUSION EQUATIONS 

The system of interest is defined as a quasi-two- 
dimensional ferromagnetic metal layer rolled out in the 
xy-plane. Two asymmetric interfaces provide a confine- 
ment in z-direction, along which the potential gradi- 
ent generates a Rashba spin-orbit coupling. Therefore 
a single-particle Hamiltonian for an electron of momen- 
tum fc is (fi = 1 is assumed throughout) 

H ^ \- aa- ■ (k X z) + -A^cO- m + H' (1) 

2m 2 

where & is the Pauli matrix, m the effective mass, 
and m the magnetization direction. The ferromag- 
netic exchange splitting is given by A^c and a repre- 
sents the Rashba constant (parameter). Hamiltonian 
— J2f=i V{f ^ fj) sums the contribution of the non- 
magnetic impurity scattering potential V{r) localized at 

To derive a diffusion equation for the non- 
equilibrium charge and spin densities, we employ Keldysh 
formalismi^ Using Dyson equation, in a 2 x 2 spin space, 
we obtain a kinetic equation that assembles the retarded 
(advanced) Green's function (G^), the Keldysh com- 
ponent of the Green function , and the self-energy 
E^, i.e., 



(2) 



where all Green's functions are full functions with inter- 
actions taken care of by the self-energies '^^'^^^ , The 
retarded (advanced) Green's function in momentum and 
energy space is 



G«(^)(fc,e) 



e-efc-(5--6(A;)-S«(^)(fc,e) 



(3) 



where = k'^ /{2m) is the single-particle energy. We 
have introduced a fc-dependent effective field b{k) — 
Axcm/2 + a{k x z) of the magnitude bk — \Axc'm/'2 + 
a{k X z)\ and the direction b = b{k)/bk- 

Neglecting localization effect and electron-electron in- 
teractions, we assume a short-range (5-function type im- 
purity scattering potential. At a low concentration and 
a weak coupling to electrons, the second-order Born ap- 
proximation is justified,— i.e., the self-energy ia^^ 



where the momentum relaxation time reads 



- « 27r 
r 



d^k' 
(2^ 



|^(fc-fc')p5(efc-efe'), 



(4) 



(5) 



where V{k) is the Fourier transform of the scattering 
potential and the magnitude of momentum k and k' is 
evaluated at Fermi vector kp. 

The quasi-classical distribution function g = 
5fc_e(T, i?), defined as the Wigner transform of the 
Keldysh function G^ {r,t;r' ,t'), is obtained by inte- 
grating out the relative spatial-temporal coordinates 
while retaining the center-of-mass ones R — {r + r')/2 
and T = {t + t')/2. As long as the spatial profile 
of the quasi-classical distribution function is smooth 
at the scale of Fermi wave length, we may apply the 
gradient expansion technique on Eq.Q;^ which gives 
us a transport equation associated with macroscopic 
quantities. The left-hand side of the kinetic equation in 
gradient expansion becomes 



[G"']-^G^ -G'^iG'^]-^ 
^[g,a-bik)] 



+ i 



dg 
dT 



+ a{z X a),VRg\ , (6) 



where {■,■} denotes the anti-commutator. The relax- 
ation time approximation indulges the right-hand side 
of Eq.© as 



p{e, T, R)G^{k, e) - G"ik, e)p{e, T, R) 



(7) 



where we have introduced the density matrix by integrat- 
ing out the momentum k in g, i.e.. 



piE,T,R) 



1 



£k' 



2'kNo J (27r) 



■9k'AT,R)- 



(8) 



For the convenience of discussion, time variable is 
changed from T to t. At this stage, we have a kinetic 
equation depending on p as well as on g 

i[6--b{k),g] + -g+^ + l\-+a{z x &),VRg 
T ot 2 \ m 



G"ik,e)p{e)-p{e)G^ik,e) 



(9) 



A Fourier transformation on temporal variable to the fre- 
quency domain uj leads to 



ng - bk[Uk,g] = iK, 



(10) 



where fl = uj + i/r and the operator Uk = o' ■ b satisfies 
UkUk — 1- The right hand side of Eq. lfTOl) is partitioned 
according to 

1 ( k 

K = --{—+ a{z X &), Vng 

2 m 



G^{k,e)p{e)-p{e)G^{k,e) 



(11) 
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The equilibrium part is denoted by 7^^°^ while the gradi- 
ent term K^^'^ is regarded as perturbation. Functions g 
and p are both in frequency domain. We solve Eq. (jlOp 
formally to find a solution to g 

. _ . i2bl-n^)K + 2blUkKUk-nbk[Uk,K] ^ 

(12) 

An iteration procedure to solve Eq. p^ has been outlined 
by Mishchenko et al., in Ref.fliJ. We follow this proce- 
dure here: According to the partition scheme on K, we 
use K^^') to obtain the zero-th order approximation as 
g^°^ = j2f[K'^^'>{p)], which replaces g in K'-^^ to generate 
a correction due to the gradient term, i.e., K'^^'> {g'^^'>). 
We further insert back to Eq.([l2l) to obtain 

a correction given by ^[K'^^^g^^'^)], then we obtain the 
first order approximation to the quasi-classical distribu- 



tion function, 

gW ^g(") +^[KW{g("))]. (13) 

The above procedure is repeated to any desired order, 
i.e., 

gin) +^[i'(l)(g("-l))]. (14) 

In this paper, the second order approximation is suffi- 
cient. The full expression of the second order approxima- 
tion is tedious thus not included in the following. A diffu- 
sion equation is derived by an angle averaging in momen- 
tum space, which allows all terms that are odd order in ki 
{i = x,y) to vanish while the combinations such as kikj 
contribute to the averaging by a factor kp6ij, given kp 
the Fermi wave vector,^- Further more, a Fourier trans- 
form from frequency domain back to the real time brings 
a diffusion type equation for the density matrix. 



■(i X (t) • p{z X &) + iC[z X (T, Vp] — B {z X a, Vp} 



F [(m X V)zp ~ o- ■ m'Vp ■ {z x &) — {z x a) ■ V pa ■ m] 

— i— (o- ■ mpa- - m- p) - iAj;c[o- • m, p] - 2R {a -m.^mx V)zp} ■ 

XC 



(15) 



r 



where all quantities are evaluated at Fermi energy ep. 
In a two-dimensional system, the diffusion constant D = 
Tv'^p/2 is given in terms of Fermi velocity vp and mo- 
mentum relaxation time r. The renormalized exchange 
splitting reads Kr^c = {/^xcl2) / {4,^'^ + 1) where = 



and 



{^iJA + a'ki)T 



C 



akpvpT 



The other parameters are 
aAxcVpkpT^ „ aA^ 



{4e + 1) 

1 2a'^kl' 



4e 



B = 



2(4^2 
2a^kl 



-1) 
+ 1 ' T,, 



R 



2{Ae 

A2 T 



1)' 



1 



Txc is the relaxation time due to the so-called D'yakonov- 
Perel mechanism.— Equation (|15p is valid in the dirty 
limit ^ <C 1, which enables the approximation H-4^^ k, 1. 
Charge density n and the non-equilibrium spin density S 
are introduced by the vector decomposition on the den- 
sity matrix p = n/2 -\- S ■ a . In a real experimental 
setupr^ii^'^S spin transport in ferromagnetic layers suf- 
fers from random magnetic scatterers, for which we intro- 
duce an isotropic spin-fiip relaxation S/r^f phenomeno- 
logically. 

Eventually, we obtain a set of diffusion equations for 
the charge and spin densities, i.e.. 



BV,-S 



+ TVz-mn + RVz-m{S -m), (16) 



at T|| " Tj_ 



'II 

AxcS X m 



1 



-m X (S X m) 



+ BV^n + 2CVz X S + 2R{m ■ V ,n)m 

+ r[m X (V^ X S) + X (m X 5)] , (17) 

where Vz = i x V. The spin density = SxX + Syij is 
relaxed at a rate 1/t|| = \/txc + ^l^sf while S'j_ = S^z 
has a rate 1/tj_ = 2/rj,c + !/'''«/• 

For a broad range of the relative strength between spin- 
orbit coupling and the exchange splitting, i.e., akp / Axe, 
Ea. ([T5)) and Eq. p7|) describe the spin dynamics in a fer- 
romagnetic layer. When the magnetism vanishes {Axe = 
0), the B-term provides a source that generates spin den- 
sity electrically;2i24 On the other hand, when the spin- 
orbit coupling is absent (a — 0), the first two lines in 
Eg. dlTp describe a diffusive motion of spin density in 
a ferromagnetic metal, which, to be shown in the next 
section, agrees excellently with early results in the cor- 
responding limit 1^ C-term describes the coherent pre- 
cession of the spin density around the effective Rashba 
field. The precession of the spin density (induced by the 
Rashba field) around the exchange field is described by 
the F-term, thus a higher order (compared to C) in the 
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dirty limit for F = IS.xctC /2. The i?-term contributes to 
the magnetization renormahzation. 



III. SPIN DIFFUSION IN A FERROMAGNET 



The last equation can be solved exactly 



5S 



1+^2 



Pf Pf,. 
— m X [je ■ V)m + — (je • v)m 

e e 



-som X 



dm 
'dt 



9m 



(21) 



Spin diffusion in a ferromagnet has been discussed ac- 
tively in the field of spintronicsi^"— In this section we 
show explicitly that, by suppressing the spin-orbit cou- 
pling, Eq. ()17p describes the spin diffusion equation in the 
corresponding limits. 

In the present model, vanishing Rashba spin-orbit cou- 
pling means a = 0, then Eq. pT)) reduces to 



= D\/'^S+—mx S 
at TA 



1 



1 



where ta = 1/ ^xc is the time scale of the coherent pre- 
cession of the spin density around the magnetization. 
This equation differs from the result of Zhang et al.j^ 
only by a dephasing of the transverse component of the 
spin density that is set by the time scale T^c- In a fer- 
romagnetic metal, we may divide the spin density into 
a longitudinal component that follows the magnetization 
direction adiabatically, and a deviation that is perpendic- 
ular to the magnetization, i.e., S — SQim -\- SS where sq 
is the local equilibrium spin density. Such a partition, 
after restoring the electric field by V — ^ V + eEd^ , gives 
rise to 



where <7 = ta(1/ts/ + l/Txc) and the electric current 
jig = e'^riTE/m is given in terms of electron density 
n. Apart from the inclusion of the dephasing of trans- 
verse component as implemented in parameter the non- 
equilibrium spin density Eq. (PT|) agrees excellently with 
Eq.(8) in Ref.flij. 

Given the knowledge of the spin density, the spin 
torque, defined as 



m X (5 X m), (18) is given by 



1 



1+^2 



—m x6S + T^SS, 

TA -Lxc 



dm dm 

vso^rr + psom X 



(22) 



at 



dt 



V — {je ■ V)m - P — m X (je • V)m 

e e 



(23) 



where rj — 1 + ^ta/T^c and /3 = ta/t^/. Assum- 
ing a long dephasing time of the transverse component 
(i.e., Txc cxd), then « 1 and Eq. P5)) reproduces the 
Eq. (9) in Ref. [2g] . On the other hand, a short dephasing 
time (of the transverse component) enhances parameter 
?7 therefore increases the temporal spin torque (i.e., the 
first term in Eq.(| 



9 ,c 9 
w;SS + ^som 
at dt 

= soDV^m + DW^SS + DePpAfpE ■ Vm 
SS Som 5S 

xc 



IV. ELECTRICALLY GENERATED SPIN 
DENSITY 



where the magnetic order parameter is allowed to be spa- 
tial dependent, i.e., m — m{r,t). The energy derivative 
is treated as d^S « PpN'F'm given Pp the spin polariza- 
tion and Afp the density of state, both at Fermi energy. 

For a smooth magnetic texture in which the character- 
istic length scale of the magnetic profile is much larger 
than the length scale for electron transport, we discard 
the contribution DV^SS.^^ The diffusion of the equilib- 
rium spin density follows soDV^m w sam/rsf. In this 
paper, we retain only terms that are first order in tem- 
poral derivative, which simplifies Eq. p^ to 



(20) 



The effect of an electrically generated non-equilibrium 
xcm X OiS, [ y) gpjj^ density due to spin-orbit coupling^ can be ex- 
tracted from Eq. (jl7p by setting exchange interaction zero 
(i.e.. Axe = 0). Retaining D'yakonov-Perel as the only 
spin relaxation mechanism and letting Tsf = oo, Eq.(|17|) 
ends up in 



_ J_m X 6S+ (— + T^]6S = 



TA 



-SQ—m + DePpNpE ■ Vm. 
dt 



DV^S {S + S,z) 

+ 2C{z xV) X S + B{z xV)n = (24) 

which reduces to the results in the well-known spin Hall 
effect.— In the case of an infinite medium along 
transport direction, i.e., ^-direction, Eq.([24l) gives rise 
to a solution to the spin density 

S =TxcBeE—ny = ^^y, 

ep TTVp 

where only the linear term in electric field has been re- 
tained. On the right hand side, we have used the charge 
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density in a 2D system n = fc|,/(27r) and introduced the 
parameter C. = akpT as used in Ref. ^24] - 

In the following sections, we explore the spin torque in 
the presence of both exchange and Rashba field in an in- 
finite medium. The primary focus is on two cases: Weak 
and a strong spin-orbit coupling, when comparing to the 
magnitude of exchange splitting. In general, Ea. (|17l) is 
applicable through a broad range of relative strength be- 
tween spin-orbit coupling and exchange splitting. A full 
scale numerical simulation on the diffusion equation is 
beyond the scope of this paper, we refer the readers to 
Ref. 1231 for further interests. 



V. WEAK SPIN-ORBIT COUPLING 

A weak Rashba spin-orbit coupling implies a small 
D'yakonov-Perel relaxation rate I/t^c a^, such that 
Txc ^ Tsf , TA , which allows spin relaxation to be dom- 
inated by random magnetic impurities. In this regime, 
when comparing to the magnitudes of C and F, the con- 
tribution from B and R are at a higher order in a, thus 
to be disregarded. We consider a stationary state where 
dS/dt = 0. An electric field applied along i;-direction, 
i.e., E = Ex. In an infinite medium,^" all the spatial 
derivatives vanishes (V — > 0) and the dynamic equation 
reads 

— y. S + X (S' X m) + -!— 5 

T"A Tsf 

=2eECy x d^S 

+ eET[yx{mxdeS) + mx{yxd,S)]. (25) 

In addition to the spin density induced by exchange 
splitting, a weak spin-orbit interaction leads to a devi- 
ation in spin density that can be considered as a per- 
turbation. Therefore, we may well apply the partition 
S = S± + S\im to separate the longitudinal and the 
transverse components. Ea. (p5| is thus reduced to 



-m X Si 



TA 



1 ^ 1 

77^*i 

Tl_ Tsf 

= -2eECPFAfFy X m 

-eErPpMpm x{yxm) (26) 

where l/Ti = l/Tj-c + I/t^/ and we have again em- 
ployed the approximation on the energy derivative d^S ~ 
PpMpTn and replaced the energy derivative of the 
charge density by the density of states at Fermi energy 
(i.e., d^n ~ n/ep = AV). We solve Eq. ([26l) to obtain a 
solution to the non-equilibrium spin density 



5, = 



TA 



eEPpMp [(2C + <rF)m x{yxm) 



1+^2 

{T -2<;C){yxm)] 



(27) 



and S\\ =0. In Ea. (l?f|) . the second component, oriented 
along the direction y x m, is actually perpendicular to 



the plane spanned by the magnetization direction and the 
effective Rashba field (along y), which, as to be shown 
below, contributes to a Rashba torque that fulfils the 
symmetry described in a recent experiment. The defi- 
nition Ea.(l22]) leads to a general expression for the spin 
torque 



T —T±y X m + T\\'m x (y x m). 



(28) 



which consists of an out-of-plane and an in-plane com- 
ponents with magnitudes determined by 



T.^^|^(2.C + /^F), 



-II 



1+^2 

eEPpMp 



(77F-2/3C). 



(29) 
(30) 



The plane is defined by the magnetization direction m 
and the direction of the effective Rashba field that in the 
present setting is aligned along y-direction. Note that 
the sign of the in-plane torque, Eq. ([30]), can change 
depending on the interplay between spin relaxation and 
precession. 

To compare directly with the results in Rcf. 10], we al- 
low the spin relaxation time Tsf — > oo, therefore /3 « 0. 
We also consider the transverse dephasing time to be 
infinite 3^1^ Under these assumptions, rj k, \ and <j w 
and we have Tj_ « 2eEPFNFC and T\\ w eEPpMpT. In 
the dirty limit, F ^ C due to A^cT <^ 1, therefore mak- 
ing use of the relation for the polarization Pp = Axc/(-f 
and the Drude relation = c^utE /ni, we obtain an 
out-of-plane torque 



am Axe 
T = 2 jey X m. 



(31) 



which agrees excellently with the spin torque in an in- 
finite system in the corresponding limit as derived in 
Ref.JQl. 



VI. STRONG SPIN-ORBIT COUPLING 

The opposite limit to Sec|V] is a strong spin-orbit 
coupling. In this case, we consider the scenario that 
akp 3> Axe and the D'yakonov-Perel relaxation mech- 
anism is dominating, i.e., l/T^e S> l/T^f^ due to the 
fact 1/Txe oc . Therefore, it is not physical to sim- 
ply assume that the direction of spin density is domi- 
nantly aligned along the magnetization direction, as what 
is treated in the case of a weak spin-orbit coupling. A 
self-consistent solution from Eg. dlTp to the spin density 
is more justified. 

Again, as in Sec|Vl we consider an infinite system 
where an electric field E is applied at x-direction. The 
magnetization direction is left arbitrary. We approxi- 
mate the energy derivative by 9^ « I/cf- The above 
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assumptions simplify Ea. (jl7p to 



A. Validity of Rashba model 



— S X m + m X {S X m) — "^^'^^ y x S 

1 eE 
+ — {S + S'^i) = — nBy, (32) 

where a strong spin-orbit coupling renders T and i? terms 
negligible. By considering T^c ^ ta,Txc, Eq. (15^ re- 
duces to 



— S xrh+ —{S + Szz) 

TA 

2eEC . ^ 

y X 5 



fF 



(33) 



which is a set of linear equations for the non-equilibrium 
spin density. We are interested in the linear response 
regime, which implies that at the distance as defined by 
the Fermi wave length 1/kp, we have eE/kp <^ akp. 
Therefore up to the first order in exchange splitting, we 
extract the spin density from the above equation to be 



eE 

7 

CF 



X 



(34) 



where x = Txc/ta we have used the identity y x m 
nizX — TTixZ. This yields a spin torque 

amAxc , ^ 

T = je [yxm 

eep 



+Xfn X {y X m) — —rrixZ x mj . (35) 

This torque is slightly different from the weak Rashba 
limit and has a strong implication in terms of magneti- 
zation dynamics. The torque is dominated by a field-like 
torque along y, similarly to the weak Rashba case. First, 
in contrast to the weak Rashba case [see Eq. (PO)) ]. the 
sign of the in-plane torque remains positive. Secondly, 
the anisotropic spin relaxation coming from D'yakonov- 
Perel mechanism yields an additional component of spin 
accumulation that is oriented along z. The implication 
of this torque on the current-driven magnetization dy- 
namics is discussed in the next section. 



VII. DISCUSSION 

Current-induced magnetization dynamics in a sin- 
gle ferromagnetic layer has been observed in vari- 
ous structures that involve interfaces between transi- 
tion metal ferromagnets, heavy metals and/or metal- 
oxide insulators. Existing experimental systems are 
Pt/Co/AlO:,^^i^i20 Ta/CoFeB/MgOfii Pt/NiFe and 
Pt/Co bilayers.— Besides the structural complexity in 
such systems, an unclear form of spin-orbit coupling in 
the bulk and interfaces places a challenge to understand 
the nature of the torque. 



The celebrated Rashba-type effective interfacial spin- 
orbit Hamiltonian was pioneered by E. I. Rashba to 
model the influence of asymmetric interfaces in semicon- 
ducting 2DEG:^ A sharp potential drop, emerging at the 
interface (say, in the xy-plane) between two materials, 
gives rise to a potential gradient VV^ that is normal to the 
interface, i.e., Vl^ ~ 'C('")^- case a rotational symme- 
try exists in the interface plane, a spherical Fermi surface 
assumption allows the spin-orbit interaction Hamiltonian 
to have the form = a&-{pxz), where a « (^)/4m^c^. 
As a matter of fact, in semiconducting interfaces where 
the transport is described by a limited number of bands 
around a high symmetry point, the Rashba form can be 
recovered through k ■ p theoryj^ 

As far as metallic interfaces are concerned, a spin-orbit 
splitting of the Rashba-type in the conduction band has 
been observed at Au surfaces Gd/GdO interfaces ji^i Bi 
surfaces and compounds)^^ and metallic quantum wellsi^ 
The presence of a Rashba interaction in graphene^Z and 
at oxide hetero-interfaces^ has also been reported re- 
cently. It is quite interesting to notice that the sym- 
metry breaking-induced spin splitting of the conduction 
band seems rather general and might not be restricted to 
heavy metal interfaces^. 

In the case of transition metals, however, the free elec- 
tron approximation fails to characterize the band struc- 
ture accurately due to both a large number of band cross- 
ing at the Fermi energy and a strong hybridization among 
s, p and d orbitals. Density functional theory (DFT) is a 
successful tool to investigate the nature of spin-orbit in- 
teraction at metallic surfaces. For example, in Refs.[39|. 
the authors observe a band splitting that possesses simi- 
lar properties as Rashba spin-orbit interaction and decays 
exponentially away from the surf ace. '^^ Alternatively, the 
spin-orbit interaction at metallic surfaces has been ad- 
dressed using tight-binding models for the p orbitals4Si^ 
At such sharp interfaces, the magnitude of the orbital 
angular momentum (0AM) is considered to play a dom- 
inant role at the onset of a Rashba-type spin splitting. 

This finding is consistent with the long stand- 
ing work on interfacial magnetic anisotropy at 
a ferromagnet /heavy metal^ and more recently 
ferromagnetic/metal-oxide interface4^ In such systems, 
a perpendicular magnetic anisotropy arises from the 
orbital overlap between the 3d states of the ferromagnets 
and the spin-orbit coupled states of the normal metal. 
The observation of perpendicular magnetic anisotropy 
at Co/metal-oxide interfaces tends to support the major 
role of large interfacial 0AM in the onset of interfacial 
spin-orbit effects^ii^ The presence of interfacial Rashba 
spin-orbit coupling has also been shown to produce 
interfacial perpendicular magnetic anisotropy4^ 

All these previous theoretical and experimental studies 
strongly suggest that the interfacial spin splitting exists 
in the presence of a large 0AM and potential gradient. 
However, a microscopic description of realistic interfaces 
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is still missing. Although the Rashba spin-orbit interac- 
tion is a convenient Hamiltonian to extract qualitative 
behaviors, its applicability to realistic metallic interfaces 
with complex band structures remains to be tested. 



B. Spin Hall effect versus Rashba torque 

Recently, Liu et al.,^^ proposed to manipulate the mag- 
netization of a Pt/Co or Pt/NiFe bilayer using the spin 
current generated by spin Hall effect in the underlying Pt 
layer. When injecting a charge current je into a normal 
metal accommodating a strong spin-orbit coupling, the 
asymmetric spin scattering induces a transverse pure spin 
current that has the form J7 = {an /e)je x p,^ p,, where 
an is the spin Hall angle and p, is the spin direction. 
When impinging on the ferromagnetic layer deposited on 
top of the Pt layer, the spin current transverse to the 
local magnetization is absorbed and generates a torque 
TsHE — (&ff/e)(l — I3mx)m x (p x m) (to be called 
SHE torque thereafter). Here, bn = OLHje^^B/e is the 
spin torque amplitude where the regular spin polarization 
P is replaced by the spin Hall angle an- /3 is the non- 
adiabaticity parameter proposed by Zhang and Li2S and 
it stems from the presence of spin-flip scattering in the 
system. In the configuration adopted by Liu et at, the 
charge current is injected along x and the torque is given 
by 

TsHE = aHfiB— {m X (y X m) + I3y X m) . (36) 

e 

Note that a more realistic model should account for spin 
diffusion in Co and Pt, as discussed in Ref. [i^. An 
important conclusion is that, besides the correction in 
the case of a strong Rashba coupling, both Rashba and 
SHE produce the same type of torque, see Ea. ((28| and 
Eg. ([55]) in this article. 

Nevertheless, distinctions can be made. First, in the 
absence of the corrections due to spin-flip and spin pre- 
cession, the Rashba torque reduces to the field-like term, 
y X m, whereas the SHE torque reduces to the (anti- 
) damping term mx[yxm). This assertion must be scru- 
tinizede carefully since the actual relative magnitude be- 
tween the field-like and the damping torques depends on 
the width of the magnetic wire as well as on the detailed 
spin dynamics in presence of spin-flip and precession. 
Furthermore, for such an ultra-small system the spin-flip 
scattering giving rise to the non-adiabaticity parameter 
(/3) might be significantly different from the one mea- 
sured in a more conventional thin film. 

A second important difference arises from the fact that 
the Rashba torque arises from spin-orbit fields generated 
by interfacial currents, whereas the SHE torque is due to 
the current flowing in the bulk of the Pt layer. Therefore, 
for a constant external electric field, varying the thickness 
of Pt layer shall enhance the SHE torque, while keeping 
the Rashba torque unchanged. 



The torques as a function of the Co layer thickness is 
more difficult to foresee. Although one could claim that 
Rashba spin-orbit interaction is expected to be localized 
at the interface, where the potential gradient is large, 
numerical simulations show that the Rashba-type inter- 
action survives a few monolayers^ (which is typically the 
thickness of the Co layer under consideration). In addi- 
tion, the presence of quantum well states might also mod- 
ify the nature of the spin-orbit interaction in the ultrathin 
magnetic layer in a system such as Pt/Co/AlO^:!^ 

The same is true for the SHE torque. The injection 
of spin current into a Co layer is accompanied by spin 
precession that takes place over a very short decoherence 
length. This decoherence length has been studied experi- 
mentally and theoretically in spin valves and found to be 
of the order of a few monolayers42. In the typical case of 
3 or 4 monolayer-thick ferromagnets, the SHE torque can 
not be considered as a purely interfacial phenomenon. 



C. Magnetization Dynamics 

In Pt/Co/AlOi: trilayers, Miron et al have observed 
a current-driven domain wall nucleation,— an enhanced 
current-driven domain wall velocity2£ and a current- 
driven magnetization switchingJ^ The symmetry of the 
spin torque required to explain the experimental findings 
agree well with Rashba torque proposed in Ref. [lOl On 
a similar structure. Pi et al^'^ and Suzuki et al^ also ob- 
served an effective field torque that could be interpreted 
in terms of the Rashba torque. Recently, Liu et aP^ inter- 
preted their experiments on Pt/NiFe and Pt/Co bilayers 
using SHE in the underlying Pt layer. 



1. Magnetization switching 

According to our previous discussions, both Rashba 
torque and SHE torque have a general form T = T±y x 
m + T\\m X [y X m). The first term acts like a field 
oriented along the direction transverse current direction 
whereas the second term acts like an (anti-)damping 
term, mimicking a conventional spin transfer torque that 
would arise from a polarizer pointing to y. 

As a consequence, both Rashba torque and SHE torque 
possess the appropriate symmetry to excite the magne- 
tization of a single ferromagnet and induce switching, as 
observed by Miron et al^^ and Liu et al?^ In the case of a 
large Rashba spin-orbit coupling, the torque acquires an 
additional component that acts like an effective magnetic 
field along i, vanishing as the magnetization component 
nix is zero (see Section IVT)) . which provides an additional 
torque that helps destabilize the magnetization. 
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2. Current-driven domain wall motion 

The influence of Rashba/SHE torque on a domain wall 
can be illustrated within the rigid Bloch wall approxima- 
tion. The perpendicularly magnetized Bloch wall is pa- 
rameterized by m = (cos sin 0, sin sin 0, cos 0) where 
(j) = (j){t) and e{x,t) = 2tan-i[e(^-^<=(*))/^], where 
refers to the center of the domain wall and A is defined 
as the domain wall width. To describe the dynamics of 
a Bloch wall, Landau-Lifshitz-Gilbert (LLG) equation 



dtm — —jm X Heff + acdtm x m + t 



(37) 



has to be augmented by the current induced torque t 

T — bj'Vm — fSbjm x Vm 

+bj{T±y X m + T||Tn x (y x m) 

+ TzTUxZ X m). (38) 

The torque t is written in the most general form, where 
the first two terms are the regular adiabatic and the so- 
called non-adiabatic torques; the next two terms (ry and 
T^) emerge from the presence of Rashba and/or spin Hall 
effect and the last term appears only in large Rashba 
limit (see Sec. IVip . The magnitude of the adiabatic 
torque is bj = iisPje/e- The effective field is given by 



ff 



2A. 



■'V m + HkitIxX + H±mzZ. 



(39) 



Parameter 7 in LLG is the gyromagnetic ratio, ac is the 
Gilbert damping, A is the exchange constant, Ms is the 
saturation magnetization, Hk is the in-plane magnetic 
anisotropy and H± is the combination of an out-of-plane 
anisotropy and a demagnetizing field. The magnetization 
dynamics can be obtained readily from Eqs. ([57)) -([M l) by 
integrating over the magnetic volume 



(9(0 -I- ac 



dtXc 



— in - y) cos 0-/3 



A 



7— sm20- 



Att 



-Ti COS( 



(40) 

b_j_ 
A' 

(41) 



We observe that the in-plane torque T|| distorts the do- 
main wall texture, while the perpendicular torque t±_ 
drives the domain wall motion. The additional torque 
Tz, arising in the large Rashba limit, only contributes to 
the in-plane torque. Therefore, in the following, we will 
refer to the in-plane torque as t*^ — — Tz/2. Below the 
Walker breakdown (dtct) = 0), the velocity is given by 



dtXc 



Avr 
~2~ 



Tji cos< 



7 sin20 



Att 



OLQ- P + —{acT^ + T|p cos 



(42) 
hj 



"gA' 



(43) 



where the tilting angle (p is given by the competition be- 
tween the magnetic anisotropy, the non-adiabatic torque, 
and the Rashba/SHE torque. In the case of weak Rashba 
(jz — 0), assuming T|| = /3tj^ and omitting the correction 
to the spin precession, we recover the results of Ref. [1^ . 
When neglecting the in-plane torque and accounting for 
the perpendicular Rashba torque (t|* = 0), the Rashba 
torque only acts like an effective transverse field and en- 
hances the Walker breakdown limiti^ [see Eq. (03])]. 

Accounting for the in-plane component T|| arising ei- 
ther from corrections to Rashba torque or from the SHE, 
this torque appears to modify the domain wall veloc- 
ity. Therefore, depending on the strength and the sign 
of Rashba/SHE torque as well as on the resulting tilt- 
ing angle 0, it is possible to obtain a vanishing or even 
a reversed domain wall velocity, as has been shown nu- 
merically in Ref. [Ill and illustrated in Eq. (gl]). A full 
scale numerical investigation is beyond the scope of this 
article, but it will help understand the profound effect of 
Rashba and SHE torque on the domain wall structures. 



VIII. CONCLUSION 



Using Keldysh technique, in the presence of both mag- 
netism and a Rashba spin-orbit coupling, we derive a 
spin diffusion equation that provides a coherent descrip- 
tion to the diffusive spin dynamics. In particular, we 
have derived a general expression for the Rashba torque 
in the bulk of a ferromagnetic metal layer, at both weak 
and strong Rashba limits. We find that the torque is in 
general composed of two components, a field-like torque 
and the other (anti-) damping one. Being aware of the 
recent alternative interpretation on the current-induced 
magnetization switching in a single ferromagnet, we have 
discussed the difference between the Rashba and the SHE 
torques. While exploring the common features, we found 
that the magnetization dynamics driven by the Rashba 
torque presents several interesting similarities to that in- 
duced by SHE torque. Nevertheless, further investiga- 
tion involving structural modification of the system is 
expected to provide a deeper knowledge on the nature 
of the interfacial spin-orbit interaction as well as the 
current-induced magnetization switching in a single fer- 
romagnet. 
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